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1 Introduction 

The law of large numbers (LLN) and central limit theorem (CLT) are long and 
widely been known as two fundamental results in probability theory. 

Recently problems of model uncertainties in statistics, measures of risk and 
superhedging in finance motivated us to introduce, in [1] and [S] (see also [2], 
[3] and references herein), a new notion of sublinear expectation, called "G- 
expectation" , and the related "G-normal distribution" (see Def. [T0|) from which 
we were able to define G-Brownian motion as well as the corresponding stochas- 
tic calculus. The notion of G-normal distribution plays the same important rule 
in the theory of sublinear expectation as that of normal distribution in the clas- 
sic probability theory. It is then natural and interesting to ask if we have the 
corresponding LLN and CLT under a sublinear expectation and, in particular, 
if the corresponding limit distribution of the CLT is a G-normal distribution. 
This paper gives an affirmative answer. The proof of our CLT is short since we 
borrow a deep interior estimate of fully nonlinear PDE in which extended 
a profound result of [T] (see also [7]) to parabolic PDEs. The assumptions of 
our LLN and CLT can be still improved. But the phenomenon discovered plays 
the same important rule in the theory of nonlinear expectation as that of the 
classical LLN and CLT in classic probability theory. 



2 Sublinear expectations 

Let be a given set and let H he a linear space of real functions defined 
on ri such that if Xi, - ■ ■ ,Xn E Ti then (p(Xi, • • • ,X„) S Ti for each S 
GpoJi/(R") where Cpoiy{^) denotes the space of continuous functions with poly- 
nomial growth, i.e., there exists constants G and fc > 0, such that |<i5(a::)| < 
G(l-|-|x|'^). Ti \s considered as a space of "random variables". 
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Here we use CpoiyiMJ^) in our framework only for some technique reason. In 
general it can be replaced by Cb(M") the space of bounded and continuous func- 
tions, by Zipf,(M") the space of bounded and and Lipschitz continuous functions, 
or by i''(M" ) the space of Borel measurable functions. 

Definition 1 A sublinear expectation K onH is a functional H [—00.00] 
satisfying the following properties: for all X^Y £0. such that ¥,[\X\], IE[|1^|] < 
00, we have 

(a) Monotonicity: if X>Y then ¥\X] >¥.[¥]. 

(b) Sub-additivity (or self dominated property): 

E[x] -E[r] <E[x-r]. 

(c) Positive homogeneity: V\XX] = XV\X], VA > 0. 

(d) Constant translatability: K[X + c] = E[X] + c. 

For each given p > 1, we denote by Tip, the collection oi X E H such that 
mX\P] < 00. It can be checked (see and ^5]) that 

E[\X + < E[\X\P]^/P + E[\Y\P]^/P. 

We also have Hq C Hp for 1 < p < g < 00 and, if ^ + ^ = I7 then for each 
X £Hp and Y £ Hq we have X ■ Y E Hi and 

E[\X-Y\] < E[|X|P]i/PE[|y|«]i/«. 

It follows that Hp is a linear space and the sublinear expectation E[-] naturally 
induces a norm := E[|X|P]i/P on Hp. The completion of Hp under this 

norm forms a Banach space. The expectation E[-] can be extended to this 
Banach space as well. This extended E[-] still satisfies the above (a)-(d). But 
in this paper only the pre-Banach space Hp is involved. 

Proposition 2 Let X,Y E Hi be such that E[Y] = —E[—Y]. Then we have 

E[X + Y] =E[X]+E[Y]. 

In particular, ifE[Y] = E[-Y] = 0, then E[X + Y] = E[X]. 

Proof. It is simply because we have E[X + Y] < E[X] + E[Y] and 

E[X + Y]> E[X] - E[-Y] = E[X] + E[Y]. 
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3 Law of Large Numbers 

Theorem 3 (Laiu of Large Numbers) Let a sequence Xi,X2, ■ ■ ■ in Ti.2 be such 
that 

E[Xf] = a^, nX^X^+,] - E[-X,X,+,] =0, J = 1, 2, • • • , (1) 
where a € (0, oo) is a fixed number. Then the sum 

5„ = Xi + • • • + X„ (2) 

satisfies the following law of large numbers: 



lim 

n — *oo 



= Km E[| — 12] =0. 

n — >OG TL 
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Moreover, the convergence rate is dominated by 

n\-\']<-- 

n n 

Proof. By a simple calculation, we have, using Proposition [21 

= ^E[Sl_, + Xl] < ^{E[5^i] + E[X'J} 
n 



Remark 4 The above condition (QP can be easily extended to the situation 
E[{X, - m)'] = ^, - - A*)] = and E[-{X, - ^i){X,+J - fi)] = 0, 

for i, J = 1, 2, • • • . Ln this case we have 

I |2i 



lini E[| — -/iPl = 0. 

4 Central Limit Theorem 

We now consider a generalization of the notion of the distribution under E of a 
random variables. To this purpose we can make a set fl a linear space of real 
functions Ti. defined on fl as well as a sublinear expectation E[-] in exact the 
same way as fi, Ti. and E defined in Section [2l We can similarly define Hp for 
p>l. 

Definition 5 Two random variables, X Ti., under E[-] and Y ^ Ti under 
E[-\, are said to be identically distributed if, for each if e Cpo/j,(R) such that 
ifiiX) £ Til, we have (p{Y) G Tii and 

E[^(x)]-E[^(y)]. 
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Definition 6 A random variable X G Ti is said to be independent under ¥,[■] 
to Y = (Yi, • • • ,Yn) S if for each test function ip G Cpoiy{M"^^) such that 
(p{X,Y) e Til, we have ip{X,y) € Tii, for each y G M" and, with Tp{y) := 
¥i[ip{X,y)], we have 

E[^iX,Y)]=EmY)]. 

A random variable X e Ti.2 is said to be weakly independent of Y if the above 
test functions (p are taken only among, instead of CpoiyiJSJ^^^), 

'fi{x,y) = My) + My)x + My)x^ i'^ea{w"), 1 = 1,2,3. 

Remark 7 In the case of linear expectation, this notion is just the classical 
independence. Note that under sublinear expectations "X is independent to Y" 
does not implies automatically that is independent to X " . 

Remark 8 // we assume in the above law of large numbers that the sequence 
XiX2,--- is dynamically independent and identically distributed from each 
other and that E[Xi] = E[-Xi] = 0, ¥\Xl] < oo. Then LLN holds. 

We denote by lipi,{M.) the collection of all uniformly Lipschitz and bounded 
real functions on R. It is a linear space. 

Definition 9 A sequence of random variables {r]i}^i in Ti. is said to converge 
in distribution under M if for each (p g lipi,(M.), {£[93(77^)]}^-^ converges. 

Definition 10 A random variable ^ € Ti. is called G-normal distributed under 
E, if for each ip e /ipf,(R), the following function defined by 

u{t, x) := E[ip{x + Vt^)], {t, x) e [0, 00) xR 

is the unique (bounded and continuous) viscosity solution of the following parabolic 
PDE defined on [0, 00) x M; 

dtu - G{dl^u) = 0, u|t=o = (/J, (3) 

where G = Ga,wi<^) is the following sublinear function parameterized by a_ and 
a with < CT < ct; 

G(a) = - CT^a"), a G M, 

Here we denote := max{0, a} and oT := (—0:)+. 

Remark 11 A simple construction of a G-normal distributed random variable 
^ is to take Q = R, H = Gpoiy{R). The expectation E is defined by E[(y9] := 
^"^(1, 0), where u = is the unique polynomial growth and continuous viscosity 
solution of (0i with Lp G Gpoiy{R) = Til- The G-normal distributed random 
variable is ^(w) = u, to G ^l = R. 
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Our main result is: 

Theorem 12 (Central Limit Theorem) Let a sequence {Xi}'^^ in Ti.^ be iden- 
tically distributed with each others. We also assume that, each Xn+i is inde- 
pendent (or weakly independent) to {Xi, • • • , Xn) for n = 1,2, ■ ■ ■ . We assume 
furthermore that 

E[Xi] = E[-Xi] = 0, E[Xf] = w"^, -E[-Xf] = g^, 

for some fixed numbers Q < a_<'a < oo. Then, {^Oi^i converges in law to the 
G-normal distribution: for each e Hpf,(R), 

lim E[^(^)]=i[^(0]. (4) 

where ^ is G-normal distributed under K.. 

Proof. For a function ip e lipi,{R) and a small but fixed h > 0, let V he the 
unique viscosity solution of 

dtV + G{di,V) = 0, (t, x) e[0,l + h]x R, V\t=i+h = if, (5) 

We have, according to the definition of G-normal distribution 



V{t, x) = E[ip{x + y/l + h-t^]- 

Particularly, ^ 

l/(/i,0) =E[(^(0], Vil + h,x)=ip{x). (6) 

Since ([5]) is a uniformly parabolic PDE and G is a convex function, thus, by the 
interior regularity of V (see Wang [6], Theorem 4.13) we have 

||V^|lci+«/2,2+c([o,i]xR) < OO, for some a e (0, 1). 

We set (5 = i and So = 0. Then 

n-1 

V{1, VSSn) ~ V{0, 0) = J2^V{{z + 1)6, VSS,+i) - V{iS, Vss^)} 

i=0 

n-1 

= {['^((* + 1)*^' - ^('^^ Vss,+i)] + [v{is, Vss,+i) - v{is, Vss,)]} 

i=0 
n-1 . 

with, by Taylor's expansion, 

»i 



[dl^V{t5, VSS^ + -yPX.+iVS) - d^^V{id, VdS^)]f3dpd-fXf^iS. 
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Thus 



dtV{i5, VdS,)6 + ldl^Vii6, VdSi}Xf^,d + d,Vii6, VdSi)X,+iV6] - E[-Is] 



i=0 

<E[V{l,VSSn)]-V{0,0) 

n-1 

< 

1=0 



n— 1 -J 

^t^^'^' ^^'^^ + oC^C*^. V5Si)Xf^,5 + d,V{i5, V6Si)X,+,VS] + E[Is] 



Since E[da;V{i5, V5Si)Xi+iVS] = E[-d:,V{i5, VSSi)Xi+iVS] = 0, and 

n^dl^Viid, V6Si)Xl,6] = E[G{dl,V{iS, VdS^ 
We have, by applying dtV{iS, VSSi) + ^dl^V{iS, VSSi) = 0, 

n — 1 



i=0 

It then follows that 

ri-l 



-E[- J2 in < nV{l, VSSr,)] - V{0, 0) < /J]. 

i=0 i=0 

But since both dtV and S^^^V" are uniformly a-holder continuous in x and ^- 
holder continuous in t on [0, 1] x R, we then have |/J| < CS'^+"/'^[l + \Xi+i\ + 
It follows that 

E[|/J|] < CT^+°/2(i+E[|Xi|°]+E[|Xip+«]). 

Thus 

-C(i)"/2(1 + E[|Xir + < E[l/(1, VSSn)] - ViO, 0) 

< C(^)"/2(l +E[|Xi|" + 

As n — > oo, we thus have 

lim E[V{1, VSSn)] = V{0, 0). (7) 

n — ^oo 

On the other hand, we have, for each t, t' G [0, 1 + h] and x E M, 



\V{t,x)-V{t',x)\ = \E[ip{x + ^l + h- tO] - E[<^(Vl + h- t'Oll 



< \E[if{x + Vl + h- to - + Vl + h-t'0]\ 



< k^\{Wl + h-t-^/l + h-t'\)E[\^W 
<C^/\t^\, 
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where denotes the Lipschitz constant of ip. Thus \V{Q, 0) — F(0, h)\ < C\/h 
and, by (O, 

\nv{i.v~5Sr,)]-nv{V~6Sr,)]\ 

= \&[V{l,y/5Sn)]-W{^ + h,V5Sn)]\ < cVh. 
It follows form ^ and © that 

limsup|E[(^(^)] < 2Cy/h. 

Since h can be arbitrarily small we thus have 

lim Eb(^)]=i[(^(0]. 



Corollary 13 The convergence ^ holds for the case where ip is a bounded and 
uniformly continuous function. 

Proof. We can find a sequence {^k}'kLi ^*Pb(I^) such that ipk ^ p uniformly 
on R. By 

+ |E[^(0] - n^.m + mM^)] - ^^,(011. 

We can easily check that ^ holds. ■ 
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